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Abstract 
We provide a systematic study of charmless Bs — > PP, PV, VV decays {P and V denote pseudoscalar 
and vector mesons, respectively) based on an approximate six-quark operator effective Hamiltonian from 
QCD. The calculation of the relevant hard-scattering kernels is carried out, the resulting transition form 
factors are consistent with the results of QCD sum rule calculations. By taking into account important 
^S^ ' classes of power corrections involving "chirally-enhanced" terms and the vertex corrections as well as 

weak annihilation contributions with non-trivial strong phase, we present predictions for the branching 
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flj ' ratios and CP asymmetries of Bs decays into PP, PV and VV final states, and also for the corresponding 

"^ ' polarization observables in VV final states. It is found that the weak annihilation contributions with 

_ non-trivial strong phase have remarkable effects on the observables in the color-suppressed and penguin- 

^ ■ dominated decay modes. In addition, we discuss the SU(3) flavor symmetry and show that the symmetry 

vol 

-^ _ relations are generally respected. 



PACS numbers: 13.25.Hw, 12.38.Lg, 12.38. Bx, 11.30. Er 



I. INTRODUCTION 

The study of hadronic charmless bottom-meson decays can provide not only an interesting 
avenue to understand the CP violation and the flavor mixing of quark sector in the Standard 
Model (SM), but also a powerful means to probe different new physics scenarios beyond the 
SM l|,y]. In the past decade, nearly 100 charmless decays of Bu,d mesons have been observed at 
the two B factories, BaBar and Belle. The experimental program to study hadronic Bs decays 
has also started, with first measurements for the branching ratios of Bs — t- (j)(j), Bs — s- K^n" and 
Bs — )• K~^K~ made available by the CDF Collaboration [3-l5|. Remarkably, the first evidence for 
direct CP violation involving B'^ — )• K'^n^ and its CP conjugate mode has also been reported by 
the CDF Collaboration, with Acp{Bs -^ K+tt') = (39±15±8)% f^. Because of the possibilities 
of new discoveries, the Bs system will be the main focus of the forthcoming experiments at 
Fermilab, LHCb and Super-B factories. 

Theoretically, analogous to the Su^d-meson decays, the charmless two-body Bs-uieson decays 
have also been studied in great detail in the literature. For example, detailed estimates have been 
undertaken in the framework of generalized factorization [a, [Tj], QCD factorization (QCDF) \t- 



10|, perturbative QCD (pQCD) [ll|, bj and soft-collinear effective theory (SCET) \li,\l4\ 



Furthermore, various New Physics (NP) effects in Bg decays have also been considered in 
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17j. With the experimental developments at Fermilab, LHC-b and Super-B factories, more and 
more theoretical studies on Bg-ineson decays are expected in the forthcoming years. 

In this work, we shall reexamine hadronic Sg-meson decays based on an approximate six- 
quark operator effective Hamiltonian, which has been applied to B,, ^ -^ PP, PV, VV (where P 



and V denote pseudoscalar and vector mesons, respectively) decays [18|, ll9|. During the evalu- 
ation of the hadronic matrix elements of effective six-quark operators, the encountered infrared 
singularity caused by the gluon exchanging interaction is simply cured by the introduction of an 
energy scale fig, which plays the role of infrared cut-off without violating the gauge invariance. 
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as motivated by the gauge invariant loop regularization method [201422 1 . In general, the infrared 
cut-off energy scale runs with the final-state mass energy, we adopt fig ~ 500 MeV in VV final 
states and fig ~ 350 MeV in PP, PV final states, due to the fact that the former concerns higher 
mass energy than the latter. 

The calculation of strong phase from non-perturbative effects is a hard task, there exists 
no efficient approach to evaluate reliably the strong phases. From the study of hadronic Bu^d 
decays in our previous papers 
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19| , it has been shown that the effective Wilson coefficients and 
the annihilation contribution with a strong phase lead to a better prediction in our framework 
to explain the observed branching ratios and CP asymmetries. As for the effective Wilson 
coefficients associated with the operators, rather than considering only the power corrections to 
the color-suppressed tree topology 23(], we also attribute possible non-perturbative corrections 



parameterized by Vi and V2 to the corresponding two operator structures {V — A) x {V — A) 
and {V-A)x {V + A). 

It is interesting to note that, with the above simpUfied prescriptions, the method developed 
based on the six-quark effective Hamiltonian ahows us to calculate the relevant B to light meson 
transition form factors, and the resulting predictions are consistent with the results of light-cone 
QCD sum rules. To further test the feasibility of our framework, we shall extend our method 
to charmless Bg — )• PP, PV, VV decays. It is seen that the predictions for the branching ratios 
in the tree-dominated Bg decays are generally in good agreement among different theoretical 
approaches, while there might be big discrepancies in color-suppressed, penguin-dominated and 
annihilation-dominated Bg decays: in QCDF approach [23|, it favors big color-suppressed and 



penguin-dominated contributions, while in pQCD approach [24], it prefers a big annihilation 
contribution. In our approach, the results stand between the QCDF and pQCD. With this 
situation, it is expected that the future precise experimental datas will give us an unambiguous 
answer. In addition, we discuss the SU(3) flavour symmetry in the decays which have been 
experimentally observed, and show that in some other interesting decay channels the symmetry 
relations are generally respected. 

Our paper is organized as follows. In Section [III firstly, we briefly review the primary six- 
quark diagrams with the exchanges of a single W-boson and a single gluon, as well as the 
corresponding initial six-quark operators. Then we present the treatments of the singularities 
caused by the gluon exchanging interactions and the on mass-shell fermion propagator, as well 
as the vertex corrections and annihilation contributions. Section IIIII contains all the input 
parameters used in our calculations. In Section IIVI we give our numerical predictions and 
discussions for Bg — )• PP, PV, VV decays. Our conclusions are presented in the last section. 
Some details on the decay amplitudes are given in the Appendix. 

II. THEORETICAL FRAMEWORK 

A. Four-Quark Operator Effective Hamiltonian 

Let us start from the four-quark effective operators in the effective weak Hamiltonian. The 
initial four-quark operator due to weak interaction via W-boson exchange is given as follows for 
the B-meson decays: 

Oi = {q'ih)v-A{^jUj)v-A, q"" = u, c, q'^ = d, s. (1) 

The complete set of four-quark operators are obtained from QCD and QED corrections which 
contain the gluon-exchange diagrams, strong penguin diagrams and electroweak penguin dia- 
grams. The resulting effective Hamiltonian (for 6 — )• s transition) with four-quark operators is 



given as follows [25l |: 
Gp 



^ q=u,c 



C,{pi)0^f{fi) + C2{fi)0i'\f^) + j;c,(^)o,(/.) 



(9), 
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i=3 



+ ll.C. 



(2) 



where Ag = VqbV* are products of the CKM matrix elements, Ci{fi) the Wilson coefficient 
functions [25,], and Oi{fi) the four-quark operators: 



Oi = {qibi)v-Aisjqj)v-A , 



^{<l) 



O2 = isibi)v-A{qjqj)v- 
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O3 = {sihi)v-A^{qjqj)v-A , O4 = '^{qibi)v^A(sjqj)v-A , 

g' 9' 

O5 = (sj6i)v-^^(^jgj)v+^ , Oe = {sibj)v-A'^{qjqi)v+A , 

q' q 

o o 

o o 

Here the Fermi constant Gp = 1.16639 x 10^^ GeV~^, {q'q')v±A = Q'lfii^ i 75)9') and i,j are 
the color indices. The index q' in the summation of the above operators runs through u, d, s, 
c, and b. The effective Hamiltonian for the b ^ d transition can be obtained by changing s into 
d in Eqs. © and ([3]). 

B. Six-quark Diagrams and Effective Operators 

As mesons are regarded as quark and anti-quark bound states, the two-body hadronic B- 
meson decays actually involve three quark-antiquark pairs. It is then natural to consider the 
six-quark Feynman diagrams which lead to three effective quark-antiquark currents. The initial 
six-quark diagrams of weak decays contain one W-boson exchange and one gluon exchange, thus 
there are four different diagrams as the gluon exchange interaction can occur for each of four 
quarks in the W-boson exchange diagram, see Fig. [TJ 



FIG. 1: Four different six-quark diagrams with a single W-boson exchange and a single gluon exchange. 



The resulting initial effective operators contain four terms corresponding to the four diagrams, 
respectively. In a good approximation, the four quarks via W-boson exchange can be regarded 
as a local four-quark interaction at the energy scale much below the W-boson mass, while the 



two QCD vertices due to gluon exchange are at two independent space-time points, the resulting 
effective six-quark operators are in general non-local. The six-quark operators corresponding to 
the four diagrams in Fig. [1] are found to be 



(2^)4 (2^)4- V. v-oy/.- -^^-^^^k^^ie 

p'^ — m^ + It 



xfe(x2) /^7'1- 7"r"rigi(xi)) * (g-4(xi)r2<z3(xi)), 

pZ — 777,^^ + ie 



xfe(xi)rigi(xi)) * (g-4(x2) /^7'1- 7'^r"r2g3(xi)), (4) 

where fc and p correspond to the momenta of gluon and quark in their propagators, qi is usually 
set to be the heavy b quark, xi, X2 and X3 are space-time points corresponding to the three 
vertices. The color index is summed between 51,52 and 53,54. Note that all the six-quark 
operators are proportional to the QCD coupling constant ag due to gluon exchange. Thus the 
initial six-quark operator is given by summing over the above four operators: 

o(^)=x:os?. (5) 

i=i 
Unlike the classical four-quark effective operator, the six-quark operators used here are non-local 
with quark and gluon propagators inserted into them. 

Before proceeding, we would like to address that the above effective six-quark operators ob- 
tained with the gluon exchanging diagrams are similar to the four fermion interactions in the 
Nambu-Jona-Lasino model. As it is well known that the effective four quark interactions pro- 
posed in Nambu-Jona-Lasino model can well lead to the dynamical chiral symmetry breaking 
via quark condensate, which was motivated from a single gluon exchange operator. The non- 
perturbative effects were taken into account from the gluon coupling constant a^ by running it 
down to the low energy scale, and from the QCD confinement scale ^ ~ ^qcd — 300 ~ 400 
MeV which may be regarded as a dynamical gluon mass scale in the infrared region: 



Such a picture has successfuhy described the nonperturbative effects of QCD at low energy 
dynamics. We all believe that perturbative QCD cannot correctly deal with the hadronic decays, 
which is actually the main motivation in our paper to develop an alternative approach to treat 
the hadronic two-body decays within the framework of QCD. It is clear that QCDF approach 
cannot compute theoretically the hadronic amplitudes from the framework of QCD and it needs 
to have some inputs for the form- factors of hadronic matrix elements. It is well known that in 
the calculations of hadronic decay amplitudes based on the effective four-quark Hamiltonian, 
the perturbative contributions of QCD are characterized by the Wilson coefficient functions 
running from high energy scales to low energy scale around fi ~ 1.0 ~ 1.5 GeV, and the 
nonperturbative contributions are carried out by evaluating the hadronic matrix elements of the 
effective quark operators at low energy, where some nonperturbative effects are considered in 
the wave functions of hadrons. In our approach with non-local effective six-quark Hamiltonian, 
the treatment is similar. The perturbative contributions are included in the Wilson coefficient 
functions with an additional QCD coupling constant as due to the gluon exchanging diagram 
for obtaining effective six-quark operators, and the gluon-gluon interactions are partially taken 
into account in the running of Og from high energy scale to low energy scale around /i ~ 1.5 
GeV. The nonperturbative contributions are considered in the evaluation of hadronic matrix 
elements of nonlocal effective six-quark operators, where the additional nonperturbative effects 
are taken into account in the non-local effective quark operators and the QCD confinement 
scale /[/ ~ ^QCD ~ 300 — 400 MeV due to gluon exchanging diagram. The strong gluon- 
gluon interactions at low energy scale are effectively characterized by a dynamical gluon mass 
or infrared cut-off scale due to QCD confinement, which is similar to the Nambu-Jona-Lasino 
model. 

With the above considerations, the QCD factorization approach with six-quark operator ef- 
fective Hamiltonian enables us to evaluate all the hadronic matrix elements of two-body hadronic 
B-meson decays. For the hadronic matrix elements relevant to Bg — ?■ PP, PV, VV decays, we 
shall list them in the Appendix. 

C. Treatment of Singularities 

In the evaluation of hadronic matrix elements, there are two kinds of singularities. One 
singularity stems from the infrared divergence of gluon exchanging interaction, and the other 
one from the on mass-shell divergence of internal quark propagator. 

In general, a Feynman diagram will yield an imaginary part for the decay amplitudes when 
the virtual particles in the diagram become on mass-shell, and the resulting diagram can be 
considered as a genuine physical process. It is well known that, when applying the Cutkosky 



rule [26| to deal with a physical-region singularity of all propagators, the following formula holds: 
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(6) 



Lr ■■■'q_ 

which is known as the principal integration method, and the integration with the notation of 
capital letter P is the so-called principal integration. 

However, the Cutkosky rule may not directly be used to treat the singularities from the in- 
frared divergence of massless gluon propagator and also the light-quark propagators due to the 
confinement and dynamical chiral symmetry breaking of strong interactions. In fact, integration 
with those propagators is sensitive to the infrared cut-off for gluon and light-quark propagators, 
and diverge to infinity when the cut-off becomes zero. A modified integration with different 
parameters for different channels is used in QCDF framework [9|, |23|, while the transverse mo- 
mentum kT dominating in the zero momentum fraction is added to the propagator in pQCD 
framework [27|. In our approach, we prefer to introduce the cut-off energy scales for both gluon 



and light quark motivated from the symmetry-preserving loop regularization 2Cll-l22l| and in order 
to investigate the infrared cut-off dependence for the theoretical predictions: 
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(q is a light quark). 



(7) 



It is noted that, as the gauge dependent term fc^/cj^ can always be rewritten as linear combina- 
tions of the momenta pa on the external lines of the spectator quark, which are all on mass-shell 
in our case (as defined in Fig. [2]), their contributions are equal to zero once the equation of 
motion is used. Our results are therefore gauge independent. 




Pi = ^-^n- 



FIG. 2: Definition of momenta in B — > A/1M2 decay. The light-cone coordinate is adopted with 



[n,n 



D. Vertex Corrections and Annihilation Contributions 



As shown in Ref. [28], the CP-violating observables may be improved by adding vertex 
corrections. Furthermore, the vertex corrections were proposed to improve the scale dependence 
of Wilson coefficients of factorizable emission amplitudes in QCDF 29|]. Those coefficients are 



C2,i-1 



always combined as C2n-i + -j^ and 6*2™ + '^^' ^ , which, after taking into account the vertex 
corrections, are modified to 



C2„-l(Ai) + -^(/i) 



C2n-.if^) + %^(/.) + ^C^%^V2.-l(M2) , 
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C2„(/X) + %^(/i) ^ C2„(/.) + %^(^) + 



"^^'^)C^^^V^.„(M2) , 



(8) 



A^, ^"' ^"^"^ A^e '"' Att ' N, 

with n = 1, ..., 5, and M2 being the meson emitted from the weak vertex. In the naive dimensional 
regulation (NDR) scheme, Vi{M) are given by [9l. l29t| 

' Vi{M) = 12ln{^) - 18 + f^ dx M^) 9ix) , for i = 1 - 4,9, 10 , 

Vi{M) = I V2{M) = -12ln{^) + 6-J^dxM'^-x)g{l-x), fori = 5, 7, (9) 

^ V3(M) = -6 + /J dx (pbix) h{x) , for z = 6, 8 , 

where (/>a(x) and </)fe(x) denote the leading- twist and twist-3 distribution amplitudes for a pseu- 
doscalar or a longitudinally polarized vector meson, respectively. While for a transversely po- 
larized vector final state, (t)a{x) = 4'±{x,i/}^ and (j)b{x) = 0. The functions g{x) and h{x) used in 
the integration are given respectively as 
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2Li2(x)-ln2x-F^^- {3 + 2iTT)lnx 
1 — X 



(x -n- 1 — x) 



h{x) = 2 Li2(x) — In^ X — (1 + 2i tt) In X — (x O 1 — x) 



(10) 



To further improve our predictions, we shall examine an interesting case that vertices receive 
additional large non-perturbative contributions, namely the Wilson coefficients Oj = Cj + -j^ 
are modified to the following effective ones: 
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C.i,) + %i(/.) + ^C.%M ( v,iM2) + I4(M2; 
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aif^) + %i(/^) + ^c^%M ( V2{M2) + F2(M2; 
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A" = C.(^) + %(/^) + ^C^%M ( y3(M2) + F2(M2; 



Nr ^"' ■ An ' Nr 



(i = 1-4,9,10), 

(i = 5,7), (11) 

(i = 6,8). 



The corrections Vi{M2) and V2{M2) depend on whether the meson M2 is a pseudoscalar or a 
vector meson, and could be caused from the higher order QCD corrections and some non-local 
effects at low energy scale as shown in Fig. [31 It can be argued from a naive dimensional analysis 
that a large number of type III diagrams in Fig. 3c are in general no longer suppressed at low 
energy scale and may lead to a significant contribution at low energy scale. While a complete 
calculation of their effects at low energy scale is not an easy task and beyond the purpose of our 
present paper, we may first treat them as input parameters and will make a detailed investigation 
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FIG. 3: The diagrams in (a), (b) and (c) are loop contributions only to the effective weak vertex (type 
I), only to the gluon vertex (type II), and for both weak and strong vertices (type III), respectively. 



elsewhere. In comparing to the vertex corrections Vi{M) (i=l,2,3), there could also be in general 
three type of vertex corrections Vi{M2). While in our computation, we have only introduced two 
kinds of vertex corrections 14 (M2) and V2(M2). This is just for the simplicity of considerations 
by assuming that the additional two kinds of vertex corrections with strong phases correspond 
to two kinds of operator structures with current-current interactions {V — A) (^ {V — A) and 
{V — A) {V + A). As it was shown in our previous work[l9| that adopting Vi{P) = 26e~3\ 
V2(P) = -26, Vi{V) = 15et*, and VziV) = -15et*, both the branching ratios and the CP 
asymmetries of most B — )■ PP, PV, VV decay modes are improved, we shall take the same input 
for the Bs -^ PP, PV, VV decays. 

As for the weak annihilation contributions, most of them are from factorizable annihilation 



diagrams with the {S — P) x {S + P) effective four-quark vertex: 



i^p ' (M) ~ / dxdy 



(x(l - y)m\ - f^j + ie)((l - y)m| - m^ + ie) ' 

A^/HM)^^ [dxdy- [f^P^-^i^^-l)^vMl-y) ^^ 

J [x{l - y)m^^ - fi^g + ie}{{l - y)m% - m^ + le) 

4if^(M)~ [dxdy ^ , ^"'i^"y"^^-7^^-\^V^\ ,, 

J (x(l -y)m| -//^ + ie)((l -y)m^ -m^ + ie) 

^^ J (x(l -y)m| -/i2 + ie)((l -y)m| -771.2 + ie) 

Since the contributions of these amphtudes are dominated by the area x ~ or y ~ 1, A^^p ^(M) 
and Ag^p ^ (M) have the same sign, while Ag^p ^ (M) and Ag^p ^ (M) have a different sign from 
Ag^p ^ (M). As a result, we use the same strong phase for A^^p ^ (M) and Ag^p ^ (M), and another 
one for ^^^^'(M) and A^'^^M). 



III. THEORETICAL INPUT PARAMETERS 

The numerical predictions in our calculations depend on a set of input parameters, such as 
the Wilson coefficients, the CKM matrix elements, the hadronic parameters, and so on. Here 
we present all the relevant input parameters as follows. 

A. Light-Cone Distribution Amplitudes 

For the -B.^-meson wave function, we take the following standard form in our numerical 



calculations 



m- 



bBM = NBx'^{l-xfexp 



1 f xniB, 
'2 V ^B. 



(13) 



where the shape parameter lob^ = 0.45 GeV, and Nb is a normalization constant. 

We next specify the light-cone distribution amplitudes (LCDAs) for pseudoscalar and vector 
mesons. The general expressions of twist-2 LCDAs are 

oo 

0p(x,/i) = Qx{l-x)[l + Y,<{lACl'\2x-l)l 

n=l 
oo 

</.y(rc,/i) = 6x(l-x)[l + ^a^(/i)C3/2(2x-l)], 

n=l 

oo 

<t>l{x,^,) = 6x{l-x)[l + ^aY{f,)C'J\2x-l)], (14) 



n=l 



TABLE I: Values of Gegenbauer moments at the scale fi=l GcV taken from [31| and at /i=1.5GeV via 
a running 





M 


■K 


K 


P 


K* 





CJ 


a-L 


1.0 

1.5 


— 


0.06 ±0.03 
0.05 ±0.03 


— 


0.03 ±0.02 
0.03 ±0.02 


— 


— 


0.2 


1.0 
1.5 


0.25 ±0.15 
0.20 ±0.12 


0.25 ±0.15 
0.20 ±0.12 


0.15 ±0.07 
0.12 ±0.05 


0.11 ±0.09 
0.09 ±0.07 


0.15 ±0.07 
0.12 ±0.05 


0.18 ±0.08 
0.14 ±0.06 


af 


1.0 
1.5 


— 


— 


— 


0.04 ±0.03 
0.03 ±0.03 


~ 


— 


T 
02 


1.0 
1.5 


- 


- 


0.14 ±0.06 
0.11 ±0.05 


0.10 ±0.08 
0.08 ±0.06 


0.16 ±0.06 
0.13 ±0.05 


0.14 ±0.07 
0.11 ±0.05 



and those of twist-3 ones are 



1, (paix,^) = 6x(l - X), 



Mx,fi) = 3[2x-l + ^a^'^(/i)P„+i(2x-l)] 



n=l 



(p+ix) = 3(1 -x)^ 0_(x) = 3x2, 



(15) 



where C„(x) and Pn(x) are the Gegenbauer and Legendre polynomials, respectively. The shape 
parameters of light mesons are taken from 3l(] and listed in Table HI 

The parameters in TableUare given at the scale /i=1.0 GeV and ^=1.5 GeV, where the values 
at /u=1.0 GeV are taken from ref.[31l], which should run to the physical scale in the B meson 



decays with /u ~ ^TlKocD^h- In our numerical calculations, we take ^i = 1.5 ± 0.1 GeV which 
is corresponding to Kqcd ^ 288l^8MeV evaluated from the data as{M:,) = 0.1172 ± 0.002. It 
is noted that LCDAs of light mesons become much closer to their asymptotic forms (all shape 
parameters become zero) when the scale runs to higher values. 



B. Decay Constants and Other Input Parameters 



For decay constants of various mesons and other hadronic parameters, we list them in Table HIl 



As for the CKM matrix elements, we shall use the Wolfenstein parametrization 



A = 0.798 



+0.023 
-0.017' 



A 



2252+°-°°°S3 n 



values of four parameters [3: 
0.340 ± 0.016. 

In our numerical calculations, the running scale is taken to be 



0.141 



+0.035 
-0.021' 



with the 
and f/ = 



^ = 1.5 ± O.lGeV ~ y/2AQCDmb. 



(16) 



The scale of as{p) in the six-quark operator effective Hamiltonian is also taken at /i = 1.5 GeV. 
The mass of b quark used here is the running mass at ^u = 1.5 GeV and evaluated following the 
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TABLE II: The hadronic input parameters [35| and the decay constants taken from the QCD sum 
36l l37| and Lattice theory [38] . 



rules 



TB, 



TUB, 



nib 



nit 



md 



L472ps 


5.366GeV 


4.4GcV 


173.3GeV 


4.2MeV 


7.6MeV 


L5GeV 


rris 


m^± 


m^o 


ruK 


TOpO 


mp± 


m^j 


0.122GeV 


O.WOGeV 


0.135GeV 


0.494GeV 


0.775GcV 


0.775GeV 


L7GeV 


TO0 


m^.± 


rriK'O 


^i^ 


m 


/b. 


u 


l.SGeV 


300MeV 


0.78GeV 


L02GeV 


0.892GeV 


0.230GeV 


0.130GeV 


Ik 


fp 


/- 


Ik- 


u 


fh 




0.16GeV 


0.216GeV 


0.187GeV 


0.220GeV 


0.215GeV 


0.185GeV 


0.186GeV 


Jp 


J UJ 












0.165GeV 


O.lSlGeV 













framework in [3^] as: 



mg(/i) 


= ^(a.(/^))"i 


9 ' 


7^(aJ 


/a^\7o/^o 
V vr / 


1- 



TT 



(17) 



The definition of Ci can be found in [34]. Numerically, we find that mfe(/i) ~5.44 GeV at NLO 
when ^ = 1.5 GeV. 

In addition, the infrared cut-offs for gluon and light quarks are the basic scale to determine 
annihilation diagram contributions (the smaller jig, the larger contributions). From our previous 
work [19], it is noted that reasonable predictions can be made when taking jiq = fig=0.37 GeV. 
Here we shall use the same values for Hg and Hg. 

C. Form Factors 

As it is known that the transition form factors for the so-called factorizable contributions in 
the usual four quark effective Hamiltonian approach have to be provided from outside in QCDF 
and SCET (such as by resorting to QCD sum rules or lattice QCD). The method developed 
based on the six quark effective Hamiltonian allows us to calculate the relevant transition form 
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TABLE III: The Bs — > P, F form factors at q-^ ~ in QCD Sum Rules, Light Cone and present work, 
where the errors stem mainly from the uncertainties in the global parameters fiscaie ~ 1.5±0.1GeV, /ig = 
0.37 ± 0.037GeV, and the shape parameters of light mesons. Within their respective uncertainties, our 
predictions are consistent with the results from the other approaches. 



Mode 


F(0) 


QCDSR[39] 


LC[40] 


LC[41] 


This work 


Bs -^ K* 


V 


0.311 


0.323 


0.285 


n 997+0.064+0.003 
"■^^'-0.037-0.002 




^0 


0.360 


0.279 


0.222 


n 9on+0-082+0.013 
U.ZOU_o. 043-0. 008 




Ai 


0.233 


0.228 


0.227 


n 178+0.046+0.002 
"■-L'°-0. 027-0. 002 


Bs -^ 4) 


V 


0.434 


0.329 


0.339 


n 9r,Q+0. 080+0. 006 

U.ZdJ_Q g3g_QgQg 




Ao 


0.474 


0.279 


0.212 


n 01 1 +0.096+0.014 




Ai 


0.311 


0.232 


0.271 


n IQ/I +0.052+0.004 
U.iy'i_o. 028-0. 002 


Bs^ K 


Fo 




0.290 


0.296 


n 9«n+0-053+0.007 
U.zou_o. 031-0. 003 



factors via a simple factorization approach. They are calculated by the following formalisms 

4:1Tas{fl)CF^FMiM2 I 






>Afi 



Ncm% Fm2 



'LL 



'■iBs){Mi,M2=P), 



V 



Bs^Mi 



^TTas{lj)CF rT.FhhM2tv ^ 
AT ^"i T? LL,± y^sj 



^BsK^Bs +'mMi, 



mM2 ("1^ 



A, 



Bs-^Mi 




A 



47ra^(/x)CF ^FMiM2 

lM2f-g 



Ncm\^FM2 



(5s) (Ml 



B^ - "T-MiW-Ma) 
y,M2=P), 



(Mi,M2 = y), 



B^^Mi _ ^'^(^s{fJ')CFrT.FMi 



Ncm%^FM2 



T, 



m 



Bs 



LL,// 



mM2 {triBs + ^Mi ) 



with: 



1 



'LL,. 



m 



LL, 



Tll, 



Cf 



N} 



(Mi,M2 = y), 



1 



(18) 



2N,, 



where the amplitudes T^^ ^ '^ are given in the appendix of [l9(](see eq. (A36)). Before giving 
predictions of the observables in Bg — >■ PP, PV, VV decays, we first present our numerical results 
for the form factor at g^ = in Table IIIII For a c omp arison, we also list the results calculated 
from the QCD sum-rules, light-cone sum rules 39144 ll| . In ref.^41i], the heavy bottom quark was 
treated based on the heavy quark effective field theory (HQEFT) via l/mq expansion. Within 
their respective uncertainties, the results obtained in our present approach are consistent with 
the ones from all other approaches. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we shall classify the 24 channels of Bg decays into two light mesons according 
to the final states, and give our predictions for the branching ratios, the CP asymmetries, 
and the longitudinal polarization fractions. Since there are only a few data available for Bg 
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TABLE IV: The branching ratio (in units of 10^^) and direct CP asymmetries (in %) in Bs — ^ PP 
decays. The central values are obtained at Hq = /ig=0.37GeV, the first error stems from the uncertainties 
in the global parameters fiscaie = 1-5 ± O.lGeV, /i^ — 0.37 ± 0.037GeV, and the second from the shape 
parameters of light mesons. NLO'^-'^ and NLO'^-^-' (0°) stand for results with "NLO correction+effective 
Wilson coefficients" and "NLO correction+effective Wilson coefficients+annihilation with strong phase" , 
respectively. 



Mode 


DatafS, 42] 


This work 




NLO 


NLC^-f-f 


NLO^^^C-IO") NLO=J^J^(5°) 


NLO'=J^J^(20°) 


Bs -> w+K- 


5.0 ±1.25 


7.7 


7.0 


7.0 7.1 


7.2 


Bs -^ ttORO 


- 


0.2 


1.1 


1.1 1.1 


1.1 


Bs -^ K+K- 


24.4 ±1.4 ±3.5 


20.8 


20.5 


17.5 22.0 


26.2 


Bs -> K'^K^ 


- 


22.6 


20.7 


17.5 22.3 


26.9 


Acp(7r+K-) 


39 ±17 


24.3 


28.5 


29.3 27.8 


24.7 


AcpinORO) 


- 


77.6 


61.6 


56.0 64.4 


72.1 


Acp{K+K-) 


- 


-14.4 


-15.4 


-18.2 -14.2 


-10.8 


Acp{K°K°) 


- 















decays, we shall make comparisons with other theoretical predictions. The comparisons with 
the current experimental data, if possible, are also made. In addition, we shall discuss the 
SU(3) flavour symmetry in the decays which have been experimentally observed and also in 
some other interesting ones. According to different decay modes, we shall give our predictions 
for the observables one by one. 

A. -Bs ^ PP decays 



This type of decays has been discussed in previous paper [l^. In this work, we will reexamine 
them and shed light on the influences of effective Wilson coefficients and annihilation contri- 
bution with a strong phase. The resulting branching ratios and CP asymmetries of Bg — )• PP 



decays are listed in Table IIVI In order to have better test on our theoretica 



framework, we also 



list the most recent predictions based on QCDF with strong phase effects 23|] and the predictions 
from pQCD [1^ approach in Table |Vl The first theoretical error in our calculations is referred 
to the global parameters of running energy scale fiscaie and the infrared energy scale Hg, and the 
second one is from the shape parameters of light mesons. 



1. B, -^ K~n+,K"n' 



0.^0 



The tree-dominated Bs — )• K^tt^ decay has sizable branching ratios of order 7.6 x 10~® 
at LO in our framework, which is above the experimental result (5.0 it 1.25) x 10^®. When 
considering the NLO contribution, one can see that the contribution would further worsen the 
prediction. After including the strong phase effect which cancels out the NLO contribution, then 
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TABLE V: Comparisons of predictions between our framework and other methods in Bs — > PP decays. 



Mode 


DatafS, 42] 


QCDF[23] 


pQCD[12] 


SCETflS] 


This work 




LO 


NLO NLOCa^-f-f.e") 


Bs -s- TT+K- 


5.0 ±1.25 


c- O+0.4+0.4 
°-'^-0.8-0.5 


y (,+3.2+0.7+0.5 
'•"-2.3-0.7-0.5 


4.9 ±1.2 ±1.3 ±0.3 


7.2 


7 7 7 1+3.2+0.6 
'•' '^^^-1.8-0.2 


Bs -^ nORO 


- 


+2.5+1.2 
^- ' -0.8-0.5 


n ic+0. 05+0. 10+0. 02 
^■J^"-0. 04-0. 05-0. 01 


0.76 ±0.26 ±0.27 ±0.17 


0.2 


r, 9 1 1+0.3 + 0.6 

U.2 i.i_o. 2-0.1 


Bs -^ K+K" 


24.4 ±1.4 ±3.5 


r,c- 9+12.7+12.5 
Z0.Z_7 2-9.1 


,o 0+4.2+7.5+0.7 
i.3.D_3 2_4.i_o.2 


18.2 ±6.7 ±1.1 ±0.5 


16.6 


20.8 22.01|:^l^lo* 


Bs -^ K'^RO 


- 


rjf. 1+13.5+12.9 

^0.1_g ;^_g 4 


1 c; 0+5.0+8.3+0.0 
iO-0_3. 8-4. 7-0.0 


17.7 ±6.6 ±0.5 ±0.6 


18.2 


22.6 22.3l5;3l3^2^ 


Bs -^ TT+TT- 


0.5 ±0.5 


n OR+O.OO+O.IO 
*J-^°-0.00-0.09 


n c.7+0. 16+0. 09+0. 01 
"J-O'-0. 13-0. 10-0. 00 




0.18 


n 9^^ n 9Q+0.01+0.26 


Bs -^ ttOttO 


- 


r, iq+0.0+0.05 
'J-J^'^-0.0-0.05 


n 00+0.08+0.04+0.01 
"J-^S-0. 07-0. 05-0. 00 




0.09 


ni9 n i9+o^oi+o^i3 

U.i2 U.i^_Qoi_o(,5 


Acpi^+K-) 


39 ±17 


20.7l-°li.| 


9. 1+3.9+3.3+2.3 
^^■^-3. 6-3. 0-1. 2 


20 ±17 ±19 ±5 


21.5 


24.3 27.8l«-0t«-? 


AcpinORO) 


- 


nf. q+17.4+26.6 
oO-J_lg.2-24.3 


CO 4+1.8+7.4+2.2 
3y-"'-4.0-1.3-3.5 


-58 ± 39 ± 39 ± 13 


1.2 


77.6 64.4j:?:01«,% 


Acp{K+K-) 


- 


7 7+1.6+4.0 
'■'-1.2-5.1 


90 O+0.9+4.9+0.8 
zo.o_g 2-4.4-1.1 


-6±5±6±2 


-15.7 


-14.4 -14.2«;11J;I 


AcpiRORO) 


- 


+0.04+0.10 
"■^-0.04-0.04 





< 10 


0.0 


0.0 0.0l01° 


Acp{n+TT-) 


- 





9+0.1+1.2+0.1 
-^■^-0.4-1.2-0.1 




5.2 


4 c 4 c: + 0.4+1.5 
4. a ^■0_o.2-0.5 


AcpinO^O) 


- 





-I 9+0.1 + 1.2+0.1 
-^■^-0.4-1.2-0.1 




5.2 


45 4 5+0.4+1.5 
4. a ^■3-0.2-0.5 



the prediction is a httle smaller than the LO one while still bigger than the data. However, our 
result is in goo d agreement with the ones obtained from the other methods, such as QCDF 



and pQCD 



12l |. In addition, the contribution of penguin operators is comparable to the one of 



tree operators, and hence the interference between the two contributions is large; as a result, a big 
direct CP asymmetry is predicted in this decay mode. Moreover, the annihilation contribution 
with a strong phase has insignificant effect on the branching ratios and CP asymmetries. 

As it is well-known that the decay Bs — )• K'tt^ can be related to B^ — >• 'k~'k~^ by SU(3) 
symmetry which implies A{Bs 
two amplitudes result in 



-^ K"7r+) ss A{Bd -^ 7r~7r+) |43|, |4j]. The relation between the 



Br{Bs -^ K-7T+) « Br{Bd -^ tt-tt+), Acp{Bs -^ K-^+) ^ Acp{Bd ^ 7r-7r+). 
These relations are satisfied experimentally. From the Table IVl and our previous paper 



(19) 



191], we 



have Br{Bd -^ ■k-'k+) = 6.6lf;^, A{Bd -^ -k'-k^) 
asymmetries, respectively. So the ratios: 

BriB, -^ K-7r+ 



26.0 



Rsr 



BriBd 



vr vr^ 



1.08, Rac 



_^2 fo^ branching ratios and direct CP 
Acp{Bs -^ K-TT+ 



1.07, 



vr vr^ 



(20) 



indicates that the SU(3) symmetry relations are satisfactorily respected in our framework. 

For Bg — )• K^t:^, its branching ratio is much smaller than the ones of Bg — )■ K~'k~^, since it is a 
color-suppressed decay mode. The contributions of effective Wilson coefficients and annihilation 
with a strong phase enhance the rate by a factor of about 5.5. As for the direct CP asymmetry, 
the NLO contribution, as well as the effective Wilson coefficients and annihilation with a strong 
phase have remarkable effects on it. It is interesting to note that the prediction of branching 
ratio in our framework is roughly consistent with the one in SCET [13], while the results of 
direct CP asymmetry have opposite signs in these two methods. 
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2. B, -^ K-K+,K"K^ 

These two decay channels are penguin-dominated, and their branching ratios are of order 
23 X 10^^ in our framework after including the effective Wilson coefficients and annihilation 
contribution with a strong phase. Such large branching rates can be easily measured at future 
LHC-b and SuperB experiments. Moreover, the annihilation contributions with a strong phase 
have remarkable effects on the branching ratios in these two decay modes. In our method, the 
prediction of direct CP asymmetry in Bg — t- K^K~ decay has the same sign as the ones in other 
methods, and the one in Bg — )• K^K^ decay vanishes because of the absence of interference 
between the tree and penguin amplitudes. It has been argued that the CP asymmetr y o f the 



decay Bg — t- K^K^ is a very promising observable to look for effects of new physics 1^, |45 |. 
For example, it is shown in [IG!] that the direct CP violation of Bg — )• K^K^ , which is not more 
than 1% in the SM, can be 10 times larger in the presence of SUSY contributions while its rate 
remains unaffected. 

In addition, the decay Bg — )■ K^K^ can be related to B^ — )• ■k'^K^ by SU(3) symmetry, 
which implies that 

Br{Bg -^ K-K+) « Br{Bd ^ t^-K+), Acp{Bg ^ K-K+) « AcpiB^ -^ tt-K+), (21) 



the first relation is experimentally satisfied 24.4 it 1.4 it 3.5 ~ 19.4 it 0.6. Our predictions 19| 
for the the two ratios are: 



BrjBg ^ K-K+) ^ ^ AcpjBg -^K-K+) 

Br{Bd^7:-K+) ^ ' ^^^ Acp{Bd ^ Tr-K+) 



^i^^= oTo 1:;^^ ^1-07, i?Acp = 7 TJ :_;;^; ^i-08, (22) 



indicates that the SU(3) symmetry relations are satisfactorily respected in our framework. 

3. Bs ^ TT+TT- , 7r°7r° 

The Bg — )• TTvr decays are pure annihilation processes. The predictions of Br(Bg — )• vr+vr 
2.8 X 10^^ and Br{Bg — > vr'^vr'^) = 1.4 x 10^^ are in good agreement with the QCDF results J23l |. 



Moreover, the direct CP asymmetry are small in these decays. Being of pure annihilation 
processes, it is meaningless to discuss the phase influence in these decays, we thus do not list 
them in Table HVl 

B. Bs^ PV decays 

We now turn to discuss the observables in Bg — )• PV decays. Since the final states are a 
pseudoscalar and a vector meson, only the longitudinal polarization of the vector meson can 
contribute. In general, there are two kinds of emission diagrams as the Lorentz structure of the 
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TABLE VI: The branching ratios (in units of 10^^) and direct CP asymmetries (%) in Bs — > nK*,pK 
decays. We use different strong phase in these decays since they are characterized by the Bg ^ V and 
Bs ^ P transition form factors, respectively. The other captions are the same as Table HVl 



Mode 


This work 




NLO 


NLO'=-f-f 


NLO'=-f-f(45°) 


NLO=-''-''(60°) 


NLO'=-''-''(75°) 


Bs -^TT+K*- 


8.2 


7.3 


7.2 


7.2 


7.3 


Bs -^ 7tOK*o 


0.2 


0.3 


0.3 


0.3 


0.3 


Acp{n-K*+) 


-27.5 


-28.3 


-18.7 


-12.8 


-6.0 


Acp{7T°K*°) 


-29.3 


66.5 


15.4 


-3.8 


-22.7 




NLO 


NLO<=-f-f 


NLO-^^^C-IO") 


NLO<=-f-f{5°) 


NLO<=-^-''(20°) 


Bs -^ p+K- 


19.7 


17.5 


17.5 


17.6 


17.7 


Bs ^ p^A-o 


0.4 


0.6 


0.5 


0.6 


0.6 


Acp{p-K+) 


17.9 


18.9 


19.2 


18.5 


16.4 


Acp{p"K") 


77.0 


-29.7 


-36.8 


-25.8 


-12.9 



TABLE VIL The same as Table |Vll but for B^ 
inputs. 



KK*,pTT decays. Here we use (5°, 60°) as the default 



Mode 


This work 




NLO 


NLO=^'^' 


(5°, 60°) (- 


-10°, 60' 


) (20°, 60' 


) (5°, 45° 


) (5°, 75°) 


Bs -s- K*-K+ 


5.8 


6.6 


7.8 


7.8 


7.8 


7.7 


7.8 


Bs -s- K-K''+ 


8.1 


7.6 


8.2 


6.5 


9.7 


8.2 


8.2 


Bs -* K'OrO 


9.0 


7.9 


8.5 


6.7 


10.3 


8.5 


8.5 


Bs -^ ROr-o 


5.6 


7.5 


7.1 


7.1 


7.1 


7.4 


6.6 


Bs -^ p-TT+ 


0.04 


0.04 


0.008 


0.006 


0.014 


0.014 


0.006 


Bs -^ n^p+ 


0.04 


0.04 


0.006 


0.003 


0.011 


0.011 


0.003 


Bs ^ ^Op« 


0.04 


0.04 


0.006 


0.003 


0.012 


0.012 


0.005 


AcpiK+K-) 


54.0 


48.8 


23.2 


23.1 


23.2 


31.4 


14.0 


Acp(K*+K-) 


-32.6 


-32.2 


-29.5 


-38.2 


-22.0 


-29.5 


-29.5 


Acp{KOK*0) 























Acp(K*OkO) 























Acp{p+7r-) 


-1.9 


-1.9 


-0.6 


-0.2 


-1.1 


-1.1 


-0.2 


Acp{7T + p-) 


-1.6 


-1.6 


-0.4 


-0.3 


-0.7 


-0.7 


-0.3 


Acpi^OpO) 


-1.7 


-1.7 


-0.6 


-0.3 


-0.9 


-0.9 


-0.3 



vector meson LCDAs is different from the pseudoscalar case. If the emitted meson is a vector 
meson, the effective Wilson coefficient is the same as that of Bg — t- PP case because they are 
both characterized by the Bg ^ P transition form factors. However, when the emitted meson 
is a pseudoscalar meson, which is characterized by the Bg ^ V transition form factors, the 
effective Wilson coefficient is smaller than that of Bg — t- PP in our framework, the detail can 
be found in Eq. (|lip . The predictions for the branching ratios and direct CP asymmetries in 
Bg — )• PV decays are listed in Tables IVIIIIXI 
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TABLE VIII: The branching ratios (in units of xlO ^) of Bg — ?► PV decays. For comparison, we also 
quote the theoretical estimates of the branching ratios in the QCDF [23| and pQCD [2J| frameworks. 



Mode 


QCDF 


pQCD 


This work 




LO NLO NLO(a<=-f^,0") 


Bs -^Tv+K*- 

Bs -^ ttOR*" 

Bs -> p+K- 

Bs ^ p°KO 

Bs ^ K'-K+ 

Bs -^ K-K*+ 

Bs -^ K*°k'^ 

Bs -> K°K*''> 

Bs -^p-7r+ 

Bs -^ TT-p+ 

Bs -> nOpO 


7 O+0.4+0.5 
'■o_o.7_o.7 
n sQ+0.80+0.84 
"■Oy-0.34-0.35 
14 7+1.4+0.9 
i4^'-1.9-1.3 

-I g+2.9+1.4 
^■^-0.9-0.6 
,, O+7.0+8.1 
^^■■^-3. 5-5.1 
,(. q+3.0+4.8 
1U.CS_2 2_4.2 

10 Pi+3-4+5.1 
1U.0_2 8_4 5 

10 1+7.5+7.7 
iU.i-3.6-4.8 
r, no+O.OO+O.Ol 
"■"^-0.00-0.01 
r, no+O.OO+O.Ol 
"■"^-0.00-0.01 
r, no+0.00+0.01 
"■"^-0.00-0.01 


+2.9+0.4+0.5 

'■0-2. 2-0. 5-0. 3 

p. 07+0.02+0.04+0.01 

"•"'-0.01-0.02-0.01 

-17 0+7. 7+1. 3+1.1 

^'•"-5.6-1.6-0.9 

n ns+0. 02+0. 07+0. 01 

"•"O-0. 02-0. 03-0. 00 

+ 1.1+2.5+0.0 

^■'-0.8-1.4-0.0 

f, pi+l. 7+1. 7+0.7 

"■"-1.5-1.2-0.3 

7 0+2. 5+2. 1+0.0 

'■■^-1.7-1.3-0.0 

. q+0. 7+2. 2+0.0 

^■■^-0.7-1.4-0.0 

n r,r,+0. 05+0. 04+0. 00 

"•^^-0.05-0.06-0.01 

n 9.4+0.05+0.05+0.00 

"•^^-0. 05-0. 06-0. 01 

n 90+0.05+0.05+0.00 

"•^■^-0.05-0.06-0.01 


7.7 8.2 7.2tlitrs 
0.09 0.2 0.3l:°)-;[lS;i 
18.7 19.7 n.etlsto.l 

0.2 0.4 0.6l!]-?lJ]-i 

5.4 5.8 7.8loi^+i:5 

o.y o.l 0.2 2 <* 2 

6.5 9.0 8.5l^-fti-^ 

4.6 5.6 7Al:lil\i 
0.03 0.04 O.OllH^tS.Si 
0.03 0.04 O.OllH^troi 
0.03 0.04 O.OllHoto.M 



TABLE IX: The same as Table Em] but for the direct CP asymmetries (in %) in the Bg -> PV decays. 



Mode 


QCDF 


pQCD 


This work 




LO NLO NLO{a=.^.^,e") 


Acpi7T-K*+) 

AcpiTT^K^O) 

Acp{p-K+) 

AcpipORO) 

AcpiK+K*-) 

Acp(K*+K-) 

Acp{KOK*0) 

Acp{K*OkO) 

AcpIp+tt-) 

AcpiiT+p-) 
Acp(vrOpO) 


9f. ,+10.8+42.2 
^"•■^-10.9-36.7 
„ 7+3.5+10.1 
^^•'-2.1-11.6 
90 Q+14.6+25.0 
^».y-14.5-23.7 
9c t;+9.2+16.3 
^^^•^-8. 8-11. 3 
-ii pi+0. 5+14.0 
^^■"-0. 4-18. 8 
r, -|r, + 0.08 + 0.05 
"•^"-0.07-0.02 
r, 4Q+0.08+0.09 
U.4y_o.o7_o.l2 
11 1+0.7+13.9 
-^^■J^-O.8-15.7 

-in 9+0.8+12.7 

-'-'J-^-O.7-12.8 




-IQ p,+2.5+2.7+0.9 
^''■"-2.6-3.4-1.4 
.7 -1+7.4+35.5+2.9 
^'•-'--8. 7-29. 8-7.0 
-, 1 9+2.4+2.3+1.2 
^^■^-2.2-1.6-0.7 
70 1+6.4+16.2+2.2 
'•^■^-11.7-47.8-3.9 
cc Q+4. 4+8. 5+5.1 
30^'5-4.9-9.8-2.5 
OR 0+2.3+2.8+1.3 
•^"■"-2.3-3.5-1.2 





-, o+0. 9+2. 8+0.1 
"^^■■^-0. 4-3. 5-0. 2 
A f;+0.0+2.9+0.6 
^•"-0.6-3.5-0.3 
-1 7+0.2+2.8+0.2 
^•'-0.8-3.6-0.1 


-24.9 -27.5 -l^.Stlitt'i 

40.1 -29.3 -3.8l^;J,l^j 

17.2 17.9 18.5t^;°l|? 
-22.4 77.0 -25.8^4:1118 
52.5 54.0 23.2ti;4t2:7 
-40.0 -32.6 -29.5tlitti 

Ol°l°o 
Otgio 
-1.8 -1.9 -0.6l«;lli;i 
-1.5 -1.6 -0.4t«;llJ| 

-1.6 -1.7 -o.fit°it\i 



1. B, -^ K*-Tr+,K*°Tr' 



*o.^o 



As for Bs — )■ K*^'ir^, K*^TT^ decays, the former is tree-dominated and the later is color- 
suppressed. Here we use smaller effective Wilson coefficients since they are both characterized 
by the Bg ^ V transition form factors. Following the same consideration in our previous work 
for B meson decays [191], we take into our calculations a strong phase 6 = 60° as the default case 
for this type of decays. The annihilation contributions with a strong phase have insignificant 
effects on the branching ratios but remarkable effects on the direct CP asymmetries in these 
decays. The Bs — s- K*^tt^ mode has the branching ratio of order 7.2 x 10^^ and big direct 
CP asymmetry, and the predictions in our framework are in good agreement with the ones in 
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QCDF [23| and pQCD [l^. As for the color-suppressed Bg — )• K*^tt^ decay, our predictions are 
different from tlie ones of the other methods (pQCD and QCDF), especially for the direct CP 
asymmetry. 

2. Bs^ p+K-,p°K° 

These two decays, being characterized by the Bs ^ P transition form factors, are tree- 
dominated and color-suppressed modes, respectively. The annihilation contributions with a 
strong phase have remarkable effects on the direct CP asymmetries, especially in Bg — )• p^K^ 
decay. For Bg — t- p'^K~ decay, our predictions are consistent with the ones in QCDF and 
pQCD. But for the color-suppressed mode, the predictions are different from each other among 
the current theoretical methods. With this situation, it is expected that the future more precise 
experimental data will give us an unambiguous answer. 

Now we shed light on the Bg — t- K*K decays. The considerations of influence from annihi- 
lation contributions with strong phases are complex since there are both Bg ^ P and Bg ^ V 
transitions, and the related predictions are listed in Table IVIIi Although these four decays are 
all penguin-dominated, they have sizable branching ratios of order 8 x 10~^ due to the fact that 
the related CKM elements (V^^V^s ~ A^) are relatively large in these decays. Moreover, it is 



interesting to note that our predictions are smaller than that of QCDF [23|] while larger than 



that in pQCD 1^. As for the direct CP asymmetries, it is large for the former two modes since 
there are strong interference between penguin and tree amplitudes, which are consistent with 
the ones in the QCDF and pQCD methods. However, for the latter two decays, the direct CP 
asymmetries vanish since there is only one type of the combination of CKM matrix elements. 

As we know that the pairs related by SU(3) symmetry are B^ — t- K*^ti^ ^ Bg — )• K*^K~ , and 
Bd —7- p~K~^, Bg —7- K*~K~^. The exact symmetry implies that the amplitudes approximately 
equal in each pair, and then the branching ratios and direct CP asymmetries are approximately 
equal, namely 

BriBg ^ K*+K-) « Br{Bd -^ K*+tt-), Acp{Bg ^ K*+K-) « Acp{Bd ^ K*+7r-), 
Br{Bg ^ K*-K+) « Br{Bd ^ p-K+), Acp{Bg ^ K*-K+) « Acp{Bd ^ p-K+).{22,) 
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Our predictions for the ratios of the above observables are 

Br{Bs -^ K*+K-) Acp{Bs -^ K*+K-) 

^""^ - Br{B, ^ K*+vr-) ^ °-^^' ^^^^ " Acp{Ba ^ K*+^-) ^ °-^' 



BrjBs -^ K*-K+) ^ AcpjBs ^ K*-K+) 

Br{B,^p-K+) ^ ■^^' ""^^^ ^cp(i?d ^ p-^+) 



These three decays proceed only through annihilation contributions. In each decay mode, 
there are both Bg ^ P and Bg ^ V transitions, which have different non-perturbative correc- 
tions as can be seen from Eq. (jlip . As a result, the influences from annihilation contributions 
with strong phase do not vanish even though these decays are pure annihilation processes, and 
the related discussions are also listed in Table IVTIl The branching ratios are at the order of 10^^ 
and the direct CP asymmetries are small in these decays. 

C. Bs^VV Decays 

There are several other observables besides the branching ratios and CP asymmetries in Bg — )• 
VV decays, such as the polarization fractions and relative phases. Naive factorization without 
annihilation contribution predicts a longitudinal polarization fraction near 100% for all Bg — s- VV 
decay modes, while the polarization anomaly in B.^ (pep (the longitudinal polarization fraction 



/l is about 35%) has been observed by the CDF [46|] experiments. Motivated by the anomaly, we 



shall study in detail the polarization, branching ratios and direct CP asymmetries in Bg — )■ VV 
decays in this section. The effects of different strong phases on branching ratio, direct CP 
asymmetry and longitudinal polarization are listed in Table |Xl For there are no experimental 
data for the most of Bg — t- VV decays, the comparisons of predictions in different theoretical 
methods are especially important, which are listed in Table IXll 

i. Bs -^ p^K*^,p+K*~ 

These two decay channels are color-suppressed and tree-dominated, respectively. From Ta- 
bles |X] and IXIl it is noted that the decay mode Bg — )■ p^K*^ has sizable branching ratio of 
order 21 x 10~^, moderate direct CP asymmetry, and very big longitudinal polarization, which 



are all in good agreement with the predictions in QCDF 23| and pQCD |l2l | methods. As 
for the color-suppressed mode Bg — t- p^K*^, it has remarkable direct CP asymmetry of 56.8% 
and longitudinal polarization of 76%, which are also well consistent with the other predictions. 
Moreover, the effects of different strong phases on branching ratios, direct CP asymmetries and 
longitudinal polarizations are insignificant in these two decays. 
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TABLE X: CP-averaged branching ratios(in units of xlO ^), direct CP asynimetries(in %) and the 
polarization fractions(in %) for Bg — > VV decays. The central values are obtained with fig—0.52GeV and 
6" = 60°. 



Mode 


Exp[35,46,47J 


This work 




NLO 


NLO^^'^' 


NLO=-f-f(45°) 


m,o^i'i'{m°) 


NLO=^'^(75°) 


Bs -^ p°K*° 


< 767 


0.6 


0.8 


0.7 


0.7 


0.7 


Bs -^ p+K— 




23.6 


21.0 


20.7 


20.6 


20.6 


Bs -^ K'-K''+ 




13.4 


12.8 


11.0 


10.4 


9.8 


Bs -^ K*OK*o 


< 1681 


15.0 


13.1 


10.6 


9.8 


9.1 


Bs -^<t><t> 


24.0 ±8.9 


22.1 


18.7 


12.1 


10.0 


7.9 


Acpip''K*0) 




66.8 


56.4 


60.8 


56.8 


50.0 


Acp{p+K*-) 




-10.1 


-10.8 


-11.3 


-9.7 


-7.2 


Acp{K'-K*+) 




20.1 


17.9 


26.3 


26.4 


24.6 


AcpiK*OK*0) 



















AcpiH) 



















fL{p°K*0) 




80 


84 


79 


77 


76 


fUp+K*-) 




96 


96 


96 


95 


95 


fL{K*-K*+) 




72 


71 


54 


48 


43 


fUK*OK*0) 




76 


72 


50 


41 


32 


fLiH) 


34.8 ±4.1 ±2.1 


71 


65 


50 


42 


31 



2. B, 



K*+K*-.K*^'K 



*o iv*0 



As for these three decay modes, they are all penguin-dominated and their branching ratios 
are at the order of 10 x 10~^. As for the direct CP asymmetry, the ones of Bg — t- K*^K*^ 
and Bg — t- (j)(j) are zero, while there is a big direct CP asymmetry in Bg — t- K*^K*^ channel. 
Moreover, recent data from the CDF Collaboration favors a huge transverse polarizations, which 
is denoted by /t = 1 — /^ in Bg — t- (j)(j) decay; our prediction is consistent with the data and 
also agrees with the one in QCDF method 23|]. In addition, the predictions of the longitudinal 
polarization fraction, which are less than 50% in Bg — )• K*~^K*~ and Bg — )• K*^K*^ are similar 
to the one in (j)(f) mode. 

One important point should be addressed is that the annihilation contributions with a strong 
phase have remarkable effects on the branching ratios, direct CP asymmetries and longitudinal 
polarizations in these decays. Especially to be mentioned, the effects of effective Wilson coeffi- 
cients and strong phase decrease the branching ratio in Bg — t- (j)(j) mode a lot. As a result, the 
prediction of Br{Bg — )• (j)(f)) is much smaller than the current data and also that in the other 
theoretical methods. This discrepancy is due to the fact that we have taken a bigger gluon 
infrared cut-off in annihilation diagram. Jig = 520GeV, following what we have done in our 
previous work for the B decays [191] . This results in a smaller branching ratio in Bg — )• (j)(j) decay 
since the bigger Jig we take, the smaller contributions to the decay amplitude. 
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TABLE XI: The comparisons in theoretical methods in Bg — > VV decays. The central values are obtained 



with n'i=0.52GeV and 



60°. The first error in our predictions arises from the varying for /i^ 



1.4 ~ 1.6 GeV, the second one stems from the shape parameters of light mesons. 



Mode 


Exp [35, 46, 4r| 


QCDF [23] 


pQCD [121 


This work 




LO 


NLO 


NLO(a'=-^-^',6»") 


Bs ^ pOA'*o 


< 767 


, q+2.0+1.7 
^•■^-0.6-0.3 


n 00+0.09+0.14+0.00 
U-'^'i-0. 07-0. 09-0. 01 


0.2 


0.6 


-1 O+0.3+0.3 
l-U-0.2-0.2 


Bs ^p+K— 




91 c+1-3+0.9 
^1-0-2.8-1.5 


o„Q+8.2+1.4+1.2 
2U.y_g 2-1.4-1.1 


22.3 


23.6 


g-i 0+13.4+2.6 
^1-0-5.2-1.8 


Bs ^ K*-K*+ 




-f,+1.0+2.3 
'•O-1.0-1.8 


(, 7+1.5+3.4+0.5 
"■'-1.2-1.4-0.2 


10.3 


13.4 


104+3.0+2.7 
10.4_25-i g 


Bs -s- K*0K*° 


< 1681 


+ 1.1+1.9 
O-0-1.4-1.7 


7 Q+1. 9+3. 8+0.0 
'■°-l. 5-2. 2-0.0 


10.9 


15.0 


Q + 3.1 + 2.5 

y-O-2. 2-2.1 


Bs -^H 


24.0 ±8.9 


,f, -,+2.6+11.3 
ID. /"-z. 1-8.8 


oc- q+8. 3+16. 7+0.0 

00.0_gg_iQ 2-0.0 


18.9 


22.1 


10 0+2.9+3.1 
10.U_2.o-2.3 


Bs -> p+p' 




n (;o+0.04+0.73 
U.O»_Q Q4_o 53 


-, ,-,+0.2+0.3+0.0 
^■"-0.2-0.2-0.0 


0.56 


0.70 


70+0.01+0.05 
^^-'^-0.01-0.05 


Bs ^ p«p» 


< 320 


n 04+0.02+0.36 
'^■■^^-0.02-0.26 


o CI +0.12+0.17+0.01 
U-01-0. 11-0. 10-0. 01 


0.28 


0.35 


qc:+0.01+0.02 
'^-■'0-0.01-0.02 


Acpip^K'O) 




46+15+10 
40-17-25 


R-l 0+3.2+17.1+4.4 
01-5-4 7-22.8-2.3 


52.2 


66.8 


cfi + 1.0+3.0 
OO.»_05_2.9 


Acp{p+K*') 




11-1-1 


a 0+1.0+1.2+0.4 
"■^-1. 2-1. 7-1.1 


-10.0 


-10.1 


+3.5+1.3 

y^ '-3. 0-1. 3 


Acp(K'-K"+) 




01+1+2 
^^-2-4 


Q q+0.4+3.3+0.3 
^■■^-0.7-3.6-0.2 


16.1 


20.1 


00 4+2.4+2.1 

^"■4-2.5-2.5 


Acp{K*Ok*0) 




r, 4+0.8+0.6 











n+0+0 


Acp{4>4>) 




r, 9+0.4+0.5 
'^•^-0.3-0.2 





n+o+o 

^-0-0 






Acp{p+P') 







o 1+0.2+1.7+0.1 
^■^-0.1-1.3-0.1 


5.8 


5.0 


c 0+1.2 + 0.4 
''■'J-2.5-0.4 


Ac Pip" P'') 







+0.2+1.7+0.1 
^■^-0.1-1.3-0.1 


5.8 


5.0 


c- 0+1.2+0.4 
O^U-2.5-0.4 


fL{p°K*0) 




gQ+4+3 
^^^-5-23 


Ar c+0.4+6.9+0.6 
^^■^-0. 3-4. 3-0. 9 


73 


80 


77+2+0 

' '-1-0 


fdp+K"-) 




92+1+1 


Qo 7+0.1+0.2+0.0 
^•^■'-0.2-0.3-0.2 


96 


96 


05+1+0 


fL{K*-K*+) 




C-r, + 3+20 

■^^-5-21 


Ao 0+5.1+2.1+3.7 
'^•^■O-4.0-2.3-1.5 


66 


72 


^stttl 


fL{K'"K*0) 




+4+22 
00_.^_25 


AQ 7+5.7+0.6+0.0 
4^^ '-4. 8-3. 8-0.0 


69 


76 


411^11 


fM<l>) 


34.8 ±4.1 ±2.1 


06+3+23 
■30-4-ig 


f-, Q+3. 6+2. 5+0.0 
Ol-y-3. 2-3. 3-0.0 


65 


71 


i2tltl 


fUp+p-) 




100 


~ 100 


~ 100 


~ 100 


~ 100 


fL{p°P°) 




100 


~ 100 


~ 100 


~ 100 


~ 100 



It is known that we can relate these three decay modes to Bd — t- K*^ p^ , K*^p^ and (pK*^ 
by SU(3) symmetry. Although our prediction of Br{Bs — )• c/x/)) = 10.0l2'i) ^^ much smaller than 
the data, it satisfies the SU(3) symmetry through the following ratio: 

BriB, - 



Rbt 



1.08. 



(25) 



Br{Bd -^ (t)K*^) 

Moreover, the predictions of huge transverse polarizations are reasonable in these Bs decays 
since the three penguin-dominated B decays also have huge transverse polarizations. As the 
data still has a big uncertainty, it is expected that more precise experimental measurements will 
be helpful to clarify such an issue. 



3. B,^p+p-,p"p° 



Now we proceed to discuss the final two decay modes, which involve only the annihilation 
contributions, thus it is meaningless to discuss the phase influence which we shall skip in these 
decays. Since there is a \/2 factor between the amplitudes of these two decays, the predictions 
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for the direct CP asymmetry, longitudinal polarization are the same but with a factor of 2 
difference in branching ratios. Our predictions for the branching ratios are smaller than those 
in pQCD jli] but consistent with QCDF \2^ method. 



V. CONCLUSIONS 

Based on the approximate six-quark operator effective Hamiltonian derived from QCD, the 
naive factorization approach has been naturally applied to evaluate the hadronic matrix elements 
for charmless two body B-meson decays. It has been shown that, when considering annihilation 
contributions and extra strong phase effects, our framework provides a simple way to evaluate 
the hadronic matrix elements of two-body Bg decays. 

For Bg — )• PP, PV final states, our predictions for the branching ratios and CP asymmetries 
are roughly consistent with those of the other theoretical methods within their respective un- 
certainties, once the effective Wilson coefficients and annihilation amplitude with small strong 
phase {O"" = 5°) are adopted. The exception here is the branching ratio of Bg — )• ■k'^K~ mode, 
which is a little bigger than the data, but it is interesting to note that the prediction is well 



consistent with the ones in pQCD method Ij]. As the current data on the branching ratio has 
large uncertainties in this mode, more precise experimental data are expected to further test our 
prediction. In the Bg — )• VV decay modes, it is noticed that there are huge transverse polar- 
ization fractions in penguin-dominated decays in our framework when considering annihilation 
contributions with a large strong phase (0° = 60°). Moreover, the prediction for the branching 
ratio in Bg — )• (f)(f) is below the current data but can be explained by the SU(3) symmetry in 
comparison with the B decays. 

Another important point should be addressed is that the method developed in this paper 
allows us to calculate the relevant transition form factors. Our predictions for Bg to light mesons 
form factors are consistent with the results of light-cone QCD sum rules and QCD sum rules. In 
this sense, we can say that our framework is reasonable from both the theoretical considerations 
and the phenomenological applications to the hadronic bottom meson decays. 

Generally, it is observed that the predictions for the branching ratios of the tree-dominated 
Bg decays are in good agreements among different theoretical methods, while there are big 
discrepancies in the color-suppressed, penguin-dominated and annihilation Bg decays. QCDF 
method 22 1 favors big color-suppressed and penguin-dominated contributions, while pQCD 
method [2J] prefers big annihilation contributions. The predictions from our method stand 
between those results given by the other two methods. It is expected that the future more 
precise experimental datas from LHC-b and Super B factories will provide a better test and 
clarify the relevant important issues. 
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Appendix: Decay amplitudes of Bg -^ PP, PV, VV modes 

As discussed in Section |TI1 the QCD factorization approach with six-quark operator effective 
Hamiltonian enables us to evaluate aU the hadronic matrix elements in two-body hadronic B- 
meson decays. Here we hst only the results for the various decay amphtudes expressed in 
terms of different topological amplitudes, which could be found in the appendix of our previous 
paper [19 1. 



The detailed calculations of the hadronic matrix elements for Bs — )• PP decays could be found 
in our previous paper [ISl . Il9l |. As for the decay amplitudes and the hadronic matrix elements 
for Bg —7- PV, VV decays, we shall list them one by one as follows. Firstly, for Bg — t- vr/j decay 
channels, we have 

AiB'g^p+TT-) = -V^gV:,E-^{Bg) + VtsV,l[2PYiBg) + lpZ'iBs)], 
A{B^,^n+p-) = -VusV:,E-^Bs) + VsV,l[2P^/{Bg)+^-PX{Bs)], 
A{B^, ^ vrV°) = l{A{B^s ^ P^^~) + ^(^s° ^ vr+p")). (A.l) 

For Bg —7- ttK* decay channels, the amplitudes are 

AiB^^^n+K*-) = VuV:,\P^'-{Bg) + '^-P^f-{Bg) + Pf-{Bg)-\p§f-{Bg)] 
-VusV:,T^'^Bg), 

A{B',^7t'K*') = -l={VuV:,[P^*-{Bs)-P§;,^{Bg)-\p^f-{Bs) + Pf-{Bg) 

-^-P§f-{Bs)] + VusV:,C^*-{Bs)}. (A.2) 

For Bg — )• pK decay channels, we have 

A{B'g^p+K-) = VtaV:,[P^^{Bg) + \pZ'\Bg) + P§''{Bg)-\p^^'\Bg)] 
-VusV:,T^P{Bg), 

A{Bl^p'K') = -^{Vt,V,l[P^^{Bg) - P§^iBg) -^P^^^iBg) + P§P{Bg) 

-^P^^P(Bg)] + VusV:,C^%Bg)}. (A.3) 
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For Bg — 7- K*K decay channels, it is given by 

+Pr'iBs) + ^P#r*(^.) - IPif^Bs) - IpEf^iBs)] 
+VusV:,[T''*''{Bs) + E'^'^'iBs)], (A.4) 

+pr\B,) - ip#r*(i3.) - Ipi^'^'iB.) - IpEf^'m], 

A{B',^K*+K-) = -VtsV,l[P^''\B,) + ^PE§''\B,) + P^^'iB,) + Pr\Bs) 
+VusV:,{T'''''{B,) + E^'^'^Bs)], 

As for the Bg — )■ VV decays, since there are three kinds of polarizations for a vector meson, 
namely, longitudinal (L), perpendicular (±) and parallel (||), the amplitudes are also character- 
ized by the polarization states of these two vector mesons. Here, we only list the longitudinal 
ones and the other ones are similar. They are given by 

A{B',^p+p-) = -VusV:,E''P{Bs) + VuV:,[2P^^\B,) + \p^^{Bs)l 

A{B', ^ p'p') = -^A{B', ^ p+p-). (A.6) 

for Bg — 7- pp decay channels, and 

A{B'g^p+K*-) = F,,y;,[P^>(s,) + ^P^|>(i?,) + p|>(i?,)-lp||V(i?,)] 

-VusV:,T^'P{Bs), 
A{B',^p'K*') = --L{y^,y;jp^V(B^)_p|>(p,)-ip^|>(P,) + p|>(i?,) 

-\p^^*'{Bs)] + VusV:,C^*P{Bs)}. (A.7) 

for Bg — 7- pK* decay channels, and 

A{B^g^K*+K*-) = -y,,y;,[P^*^*(p,) + ^pf^*^*(p,)+Pf*^*(p,) + Pf*^*(i?,) 

+Pr''\Bg) + lpi^'^\Bg) - \p^§'^\Bg) - \pE§'^\Bg)] 

+VusV:,[T^'''\Bs) + E^'^\Bs)], 

A{B^,^K*^K*^) = -yi,F;jp^*^*(i?,)-ipf|*^*(p,)+p|*^*(p,) + pf*^*(p,) 

+pr^*(B.) - lpif^*{B,) - Ipif^^B,) - '-pEr'^'mm 
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for Bg — ?• K*K* decay channels, and 

+Pf{B,) - Ipi^HBs) - Ip^t^iBs)]. (A.9) 

for Bg —7- (j)(p decay channel. 
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